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". Abstract 

o : 

An explicit realization of anyons is provided, using the three-body Calogero 

model. The fact that in the coupling domain, — 1/4 < (7 < 0, the angular 

spectrum can have a band structure, leads to the manifestation of the desired 

O ! phase in the wave function, under the exchange of the paticles. Concurrently, 

psj ; the momentum corresponding to the angular variable is quantized, exactly 

akin to the relative angular momentum quantization in two dimensional any- 

1 — 1. onic system. 

43 ■ 

I Introduction 

Fractional statistics (FS) has been widely explored in two spatial dimensions 
PP, [2], 0, where it can arise due to the Abelian nature of the rotation group. 
Furthermore, braid group replaces the permutation group in the plane leading 
to the possibility of anyons pEj-[B]- The Chern-Simons (CS) theories in 2 + 1 
dimensions describe anyonic excitations [7], [8], where an Abelian gauge field 
with CS action implements statistical interaction [TU]-[T2] . The correspond- 
ing quantum mechanical problem led to the statistical 1/r 2 interaction in the 

■^j- ■ plane, which has been extensively studied [9j. The possibility of self-adjoint 

extensions and singular but normalizable wave functions, have been expli- 
cated in models with statistical interaction |13|. In contrast to this, in one 

OO — 

dimension, statistics and interaction get non-trivially intertwined. However, 

unlike two spatial dimensions, here solvable models exist, which show the ef- 
fect of statistical interaction explicitly. The Calogero-Sutherland model with 
its 1/r 2 mutual interaction between particles [14J, has shown the realization 



X 



of Haldane statistics in one dimension in the coupling range g > [15] - pj5] . 
Here, the energy spectrum is linear with E ni — J2i n>i + y + N ^ N ~ 1 ' (a + 5), 

where n\ = 0, 1, 2.., a = \J{1 + 2g)/2 and 5 is and 1 for bosons and fermions 
respectively. Subsequently, an explicit map has been given relating two par- 
ticle fermionic wave functions to the corresponding bosonic ones [23]. Later 
on, this has been generalized to N particles in an oscillator potential [24], as 
well as for free particle scenarios. In these studies, the wave functions show 
fermionic or bosonic behavior by construction. 
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Here, we explicitly solve the three-body Calogero problem with anyonic 
boundary conditions and establish FS through the phase acquired by the 
wave functions. In the Jacobi co-ordinates, the three body Hamiltonian 
reduces to a separable equation describing the dynamics of a particle on a 
plane, in a non-central potential for which exact solutions exist. The solutions 
of the angular part are central to this analysis, as the interchange of the 
particles and their interaction are both reflected in the angular variable. 
By mapping the potential in the angular equation to the Scarf potential 
[25], which can display band spectrum for a certain range of the potential 
parameter, anyonic boundary condition can be implemented on the wave 
function. When the particle positions are exchanged, this boundary condition 
allows the wave function to pick up a phase, which is fractional for solutions 
inside the bands and integral for the band-edge solutions. We discuss the 
solutions for the radial part briefly and present the complete solutions of the 
three body Calogero model, in the parameter range — 1/4 < g < 0. Thus 
this model presents a unique situation, where the FS is established via the 
phase acquired by the wave function. In addition to this, we also show that 
in the context of anyons, the quantized momentum associated with the band 
spectrum, takes the form similar to the quantization of angular momentum 
of anyons in two dimensions [1] and the spacing of linear momentum in one 
dimensions discussed in [26] . 

In the next section, we analyze the three-body Calogero model and map 
its angular eigenvalue problem to the Scarf potential. The solutions for the 
angular and radial equations are discussed and the complete solutions for 
the Calogero model are presented. In section III, we study the effect of the 
particle interchange on the solutions. In section IV, we discuss the quanti- 
zation of the momentum in the angular variable. Here, we briefly discuss 
other three-body problems, where anyonic excitations can exist, followed by 
the conclusions in the last section. 

II Band structure problem in the Calogero model 

The Hamiltonian for the Calogero model [27] is (h = 2m = 1), 

3 Q2 Jl 

where g > —1/4. Using the Jacobi co-ordinates, 

(xi + x 2 + x 3 ) (xi - x 2 ) {xi +x 2 - 2x 3 ) 

x = 3 • X = —J5-' V = 75 ' (2) 

we factor out the center of mass of the three particles and obtain a system 
with two degrees of freedom. Next we map it on to the (r, <p) coordinates 
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using 



x = rsin , y = rcos , r 2 = \[( x i — x 2) 2 + (x 2 - x 3 ) 2 + (x 3 -£i) 2 ]. (3) 

The range of these variables is < r < oo and < <f> < 2n and the ordering 
of the particles is reflected in the variable [27J, for example the ordering 
x\ > X2 > X3 corresponds to < < it/ 3. Thus (CQ) reduces to, 



where U(r) = \ujt 2 and U(<j>) = (sin 9 f 0)2 [2E]. Using 

*{r,<f>) = ±:u(r)K{<l>), (5) 
V r 

we separate into the radial and angular equations. The solutions of the 
radial equation, 

/ d 2 ~, , (m 2 - 1/4) \ . , * 

1 +U(r) + ± -— ^ — - ] u[r) = Eu(r), (6) 



dr 2 r 2 



are 



«(r) = r m +^ exp f-^) Lf 1 (^r 2 ) (7) 
where L™(r) are the Laguerre polynomials. The eigenvalues are 

E = (21 + m + l)u, (8) 
where m 2 is the eigenvalue of the equation in given below 

In the literature, the solutions for the Calogero model have been largely 
confined to g > regime [27] -|29j. Our interest lies in the solutions of ([I]), 
when — 1/4 < g < [3U], because for g lying in this range, @ exhibits band 
structure resulting in non-normalizable solutions [23], [3"T] . 

This is expected because for values of (7 lying in this range, ( sjn 9 f ^2 is 
similar to the potential in a crystal lattice and hence, the band structure. In 
terms of the variable, we can explain this by applying the circular boundary 
condition, K(0) = K(2n). 
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As pointed out earlier, the solutions of © will play an important role in 
determining the type of statistics obeyed by the particles. Using the mapping 
30 = — , we can map ([9]) to 



7T 



d 2 



a 2 dx 2 



[sin^f) 2 j 



K(x) = \ 2 K(x), 



(10) 



where g 



9m 2 



A 2 and a is the lattice period. Owing to the 



constraint on g, we see that s should lie in the range < \s\ < ~. This band 
structure problem has been discussed in detail by Scarf in [25] . 

The most general solution for (flOl) . written in the self-matching form 
E5J,[32], is 



K(x) 



V(X 



v 



Na) ^u(x - Na) 

- + Vp- 



Uq 



(11) 



for (N — \)a < x < (N + |)a, with N = |, | Here v (x) and u(x) are the 

linearly independent solutions of (flOl) . whose explicit form, for < x < a/2, 
is as follows 



U (X) 



V[X) 



and 





( 7TX N 




2 ) 




= sin 2 








2-^1 




V a / 








( TTX S 




"2") 




= sin 2 








2-^1 




V a / 







1 |s| 

4 + y 



A 1 

2' 4 



2 2' 



A 1 |s| A 
2' 4 _ T ~ 2' 



s ; sm 



s|; sm 



TCX 

a 

TlX 

a 



(i + Vp)/(i - Vp) >p 



fka s 
i tan — 
I 2 , 



(12) 
(13) 



We point out here that for N = 1/2, the range of x corresponds to the 
particle ordering x\ > x^ > x$. The eigenvalue m 2 corresponding to (1111) is 



m 



9tt^ 



[cos 1 (sin 7r\s\ cos ka)} x 



(14) 



where the n th principal value of the inverse cosine function is to be taken. 
Here, k is the reduced wave number. For the n th band, it takes values and 
it /a for the lower and upper band-edges respectively and values in between 
them inside the band. Substituting ([TBI in (jllj) . we obtain K(x) in terms of 



Kk{x)) = exp iNka — 



ika 



( ka\ v(x — Na) ( ka\ u(x — Na) 

cos\ — — + 1 sin — - 

V 2 J v \ 2 J u 

(15) 
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where the subscript k is used to show the dependence of the solution on the 
wave vector. It can be seen easily from the above equation that K(<f>) satisfies 
the Bloch condition: 



Kk(x + a) = exp(ika)Kk(x) 



(16) 



using which, we can obtain the solutions for the entire lattice. Substituting 
k = and k = it /a in ( HUEl) we obtain the band-edge eigenvalues for the lower 
and upper band-edges as 



rrin 



9tt 2 



n 



, m ir/a 



9vr 2 



n + \s\ 



(17) 



Similarly for these values of k, K(x) in ({To']) becomes pure v (x) and pure 
u(x) for the lower and upper band-edges respectively and the hypergeometric 
functions in these solutions reduce to the Jacobi polynomials [31]. Thus for 
the n th band, the lower band-edge eigenfunction is 



ika 



Kq(x) = exp y,Nka o 



sin 



»_ Hi _i_ I 

TlX \ 2 2+4 



n 



-i cot 



TTX 



where z/i = — n + \ s\ — \ and the upper band-edge solution is 



Kz. (x) = exp iNka 



ika" 



sin 



TTX 



ki 

2 



n 



-i cot 



TCX 

a 



(19) 

where v\ = —n — \s\ — \. Applying the Bloch condition to the band-edge 
solutions (ITS]) and (fT9"l) . we see that K (x) does not acquire a phase and 
Ktl(x) changes sign when one goes from one cell to the other. We would 
like to point out here that the above solutions f TTTT) . fp~5|) . ( fTHl) and ( TT91 are 
defined in only one lattice cell. In the <fi variable this corresponds to the 
range < <fi < tt/3. Incrementing <p by <p — > <fi + pir/3, where p — 1,2, 3, 4, 5 
defines the solutions for the entire range of <p |27j and this in the x variable, 
amounts to translating from one cell to the next cell in the lattice. Thus, by 
continuously incrementing by tt/3, we can cover the whole real line in the 
x variable. 

Having described the solutions for both the radial and the angular eigen- 
value equations, we write the complete solution of (@J as 



\I>(r, 



r exp 



Lur 



-ur 



KM, 



(20) 
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where we substitute the suitable solutions -K&(0), with the corresponding m 
values, depending on the value of k being considered. For k lying inside the 
band, the complete solution for (jl]) is 

¥(r, 0) = rH exp (-^ (^r 2 ) K fe (30), (21) 

where the m values are obtained from (114p . 

Similarly for = and k = ir/a, the complete solution \l/(r, 0) is 

\ Hmo] Z^ 1 2 



¥ (r, 0) = r' m °' exp 1 L\ m ° 1 (^r 2 J K o (30) (22) 

and 

*(r, 0) = rl*"*/-l exp J 4^' (^ 2 ) K*(ty), (23) 

where m values are obtained from ( TT71) . Thus we can see that for —1/4 < g < 
0, the Calogero model has different solutions, depending on the wave vector, 
owing to the lattice like behaviour of the potential in the variable. For 
completeness, we point out here that for g > 0, the complete solutions of the 
Calogero model are a product of the Gegenbauer polynomials comming from 
the angular part and the Laguerre polynomials comming from the radial part 

III Behaviour of the wave function under particle exchange 

As mentioned earlier, the interchange of particles is reflected in the variable. 
The range < < ir/3, corresponding to the particle ordering x\ > x 2 > x 3 , 
leads to < x < a in the x variable, which corresponds to the first cell 
of the lattice. The change in particle ordering, brought about by changing 
0^0 + jwr/3 with p = 1,2,3,4,5, translates to moving from one cell to 
the next. This leads to the wave function picking up a phase as can be seen 
easily from (1151) . The type of phase acquired, namely integral or fractional, 
will decide the statistics obeyed by the system. 

For explicitness, we change the ordering from x\ > x 2 > £3 to x 2 > 
x\ > 23 which is obtained by incrementing — > + 27r/3. The range 
7r/3 < < 27r/3 corresponds to a < x < 2a which is the second cell. Thus 
an interchange of particles requires an increment in 0, which in x variable 
corresponds to a translation to the second cell. The solutions in the first and 
second cell are obtained by putting N = 1/2 and iV = 3/2 in ( fl5i) as 



K\, k(x) = cos 




n x - 2- 

v 




u(x - f) 



wo 



(24) 



6 



with 





ka\ u(x 




(25) 



Here the additional subscripts, 1 and 2, in the left hand side of ([21]) and (|2*5|) 
denote the number of the cell. Thus when we translate from the first cell to 
the second cell of the lattice, we see that the wave function picks up a phase 



which we define as the statistics parameter. O takes values lying in the range 
< < 7r for < k < it /a. It is easy to see that every time the ordering of 
the particles is changed, the wave function picks up a phase which is fractional 
when k lies in the range < k < irja. Thus for k values lying inside the 
band, the particles obey anyonic statistics. For k = 0, which corresponds 
to the lower band-edge solution (|TSI) . 6 = 0, which implies the translation 
from first cell to the second cell does not affect the wave function i.e., the 
lower band-edge solutions are symmetric implying a compliance to bosonic 
statistics. On the contrary, for k = it /a, which corresponds to the upper 
band-edge solution ( 1191) . 6 = n, showing that the interchange leads to the 
wave function being antisymmetric, which implies that fermionic statistics 
are obeyed. Moreover, substituting ka = in (THj) . we can see that the 
angular eigenvalues are functions of the statistics parameter and for a given 
band they continuously interpolate between the bosonic and fermionic eigen 
energies. This interpolation between symmetrization and antisymmetrization 
has been discussed by Wu in the context of three particles in an harmonic 
well in two dimensions |33J and by Yang et. at, in one dimension for a Bose 
gas with two-body delta-function interactions [34|. 

Thus, we can see that, the band structure in the eigenvalue spectrum 
of the differential equation in <fr leads to rich physics, where the particles 
obey three different statistics depending on the value of k. It is interesting 
to note that in [TO], the CSM (g > 0), with N-body interaction, has been 
shown to obey FS as defined by Haldane. In this case the energy spectrum 
is linear and the scale invariant energy shift is seen as the basic reason for 
the occurence of the Haldane statistics. This result holds for the three-body 
Calogero model also as its energy spectrum is linear [27] , in the g > regime. 
When we compare this to the present case, we see that for — 1/4 < g < 0, 
the spectrum given by (JHJ) in nonlinear owing to the fact that m, given by 
(fl4l) . is a non-integer, where fractional statistics manifest. The fractional 
phase is picked up by the wave function under particle exchange. We expect 



9 = ka, 



(26) 
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this result to hold when our analysis is extended to the general CSM with 
-1/4 < g < 0. 

IV Quantization of the momentum in the variable 

Exploiting the well developed theory of energy bands in solid state physics, 
we show that the momentum in the <fi variable, corresponding to the anyonic 
statistics, is quantized. This is similar to the relative angular momentum 
quantization, proposed by Wilczek, for anyons in two dimensions [1J. 

The effective crystal momentum associated with the particle in a lattice 
is defined as 

p(k) = ^ln[r(k)). (27) 

Substituting p = ztan f^) f° r r i n (O, gives r = exp(zfca), which in turn 
gives 

p = — \n(exp[ika + 2n'iri]) 
ia 

= (28) 

Here n' is a non-negative integer. Substituting ka = 6 from ( 1261) . we get 

h r 6 ,, , N 

P=-l^Z + A (29) 

a ATT 

which is the quantization condition for p in terms of the statistics parameter. 
It is interesting to note that we can extend this analysis and establish frac- 
tional statistics in other exactly solvable three-body problems discussed by 
Wolfes [35], Khare and Bhaduri [36], where the angular eigenvalue equation 
will exhibit band spectrum for a small range of potential parameters. Simi- 
lar analysis of these problems shows that inside the band, the particles obey 
anyonic statistics and at the band edges, the statistics obeyed are bosonic 
and fermionic. We briefly discuss the Wolfes potential here. 
Wolfes potential 

The eigenvalue equation for the three-body linear problem discussed by 
Wolfes m\ is 



3 d 2 3 1 

" ~pT2 + h u2 ( X i - X if + 2 9( X i - X iY 2 \ ] ^i^k(xi, X 2 , X 3 ) 

1=1 ° X i i,j=l " 

+ ( E 6 fi( X i - X i) + ( X J - X kT 3 - E ) ^i^jM X U X 2, X 2) = 0, (30) 
\i,j,k=l 
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where g > — 1/4 and / > —1/4. In the polar coordinates, this model can be 
separated to obtain 



cf 2 



dx 2 r dr 



1 + u 2 r 2 + ^-e) X (r) = 



(31) 



and 



d 2 



+ 



+ 



9/ 



A 2 K{<f>) = 0. 



(32) 



sin z 30 ' cos 2 30 

Similar to the Calogero model studied earlier, the angular equation exhibits 
band spectrum for — 1/4 < g < and — 1/4 < / < 0. Using the mapping 
30 = — and defining t = sin 2 (^), g = (s 2 — |) and / = (w 2 — wo nari 
map the angular equation (1521 to an equation similar to ( fTUl) 



4/' 



d 2 



+ 



9<? 



7T^ 



sm 



2 / 7T3; ' 



+ 



9/ 



COS 



2 ( TtX* 



- A 2 (0) = 0. 



(33) 



By substituting <j>(z) = t a (l - t) b F(t) in ([33]), it can be reduced to a hy- 
pergeometric equation, whose linearly independent solutions in terms of x 
are 



u[x) 



and 



sm 



2 /^ 



1 4. Ill 
4" 1 " 2 



1 |s| |u>| 

2 + T + T" 



cos 



f TTX\ 



\ a J 
A 1 |s 
6' 2 + ~2 



M 

2 



A 



+ ^r + ^ + 7r;l + s ;sin' 
o 



7TX 

a 



(34) 



sm 

1 

2 



7rxy 


1 1 i s 

4 T 2 








cos 2 | 


a /. 






s 


\ w \ 


A 1 


y ~ 


' ~Y ~ 


" 6' 2 



TTX 

a 



I s I 

Y 



~2~ 



A 
6 : 



s | ; sm 



7re 
a 



(35) 



Proceeding in the same way as in the Calogero case, the above two so- 
lutions can be used to write the most general solution in the self-matching 
form, for (N - 1/2) < x < (N + 1/2), given in (fTTT) . As expected, at the 
band-edges, the self- matching solution becomes pure v(x) or u(x) with the 
hypergeometric functions reducing to the Jacobi polynomials. The solution 
K(x), with v(x) or u{x) substituted from (l34"l) and ( 1351) . picks up a phase 
0, when we travel from one cell to the other. As discussed earlier, in the 
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variable, this translates to exchanging the particle positions. Not suprisingly, 
is fractional for < k < it /a and ±1 for k = and k = ir/a respectively. 
Thus, it is sufficient to conclude that the three-body problem described by 
the Wolfes potential will also show fractional statistics in the angular equa- 
tion. From above, it is clear that this analysis can also be extended to similar 
three-body interactions discussed by Khare and Bhaduri in [HS] , as the angu- 
lar equation in all these problems also displays band spectrum in some range 
of the potential parameters. Thus, the fractional statistics seems to be a 
common feature for all the three-body interactions, whose angular potential 
exhibits lattice like behaviour for some parameter range of the potential. 
VI Conclusion 

In this paper, we have analyzed the solutions of the Calogero model for 
— 1/4 < g < and we have shown that the band structure in the angular vari- 
able has interesting consequences for the statistics obeyed by the particles. 
The angular eigenvalues turn out to be functions of the statistics parameter 
6, which is the phase picked up by the wave function when two particles 
are exchanged. For 6 lying in the range < < tt, the statistics obeyed 
is anyonic. For = and = tt, the particles obey bosonic and fermionic 
statistics respectively. We have also shown that in a given band, there is 
a continuous interpolation from the bosonic state to the fermonic state. In 
addition to this, the quantization of momentum in the variable, in terms 
of the statistics parameter, is obtained. 

We have also showed that this type of statistical behaviour is common 
for the whole class of three-body problems, whose angular equation dis- 
plays a band spectrum. It will be interesting to extend this study to mod- 
els with N > 3 and finally to the iV-body CSM in the parameter regime 
—1/4 < g < 0. The possible difficulty would be the nonavailability of ex- 
act solutions in this regime. Hence, we may have to resort to the use of 
numerical techniques and this will be discussed else where. Here, we would 
like to emphasize two points. First, the anyonic phase picked up by the 
wave function is due to the exchange of the particles. Though the statistics 
here are anyon-like, these are not analogues of the braid statistics obeyed 
by anyons in two dimensions [1] and hence there are no braids associated 
with these statistics. Secondly, these statistics are entirely different from 
Haldane's fractional exclusion statistics which are based on the generalized 
Pauli's exclusion principle [2]. As mentioned earlier, such statistics is obeyed 
by the Calogero model in the g > regime and the main reason for this 
non-trivial statistics here, is the linearity of the energy spectrum. 

The intense research on Bose-Einstein condensates (BECS), has led to 
the development of experimental techniques to control and manipulate a few 
atoms in suitable traps. In light of these developments, it looks feasible for 
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the above anyon model to be realized experimentally [38J. It is interesting to 
note that there have been proposals to use anyons for quantum computation 
|37j . Hence, the realization of FS in one dimension, in the class of three body 
problems, may have useful applications in the area of quantum computation 
and in the engineering of entangled states of atoms. 
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